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ABSTRACT 

We present analysis of 4U 1626—67, a 7.7 s pulsar in a low-mass X-ray binary system, observed 
with the hard X-ray detector of the Japanese X-ray satellite Suzaku in March 2006 for a net exposure 
of ~ 88 ks. The source was detected at an average 10-60 keV flux of ~ 4 x 10 -10 erg cm -2 s" 1 . The 
phase-averaged spectrum is reproduced well by combining a negative and positive power-law times 
exponential cutoff (NPEX) model modified at ~ 37 keV by a cyclotron resonance scattering feature 
(CRSF). The phase-resolved analysis shows that the spectra at the bright phases are well fit by the 
NPEX with CRSF model. On the other hand, the spectrum in the dim phase lacks the NPEX high- 
energy cutoff component, and the CRSF can be reproduced by either an emission or an absorption 
profile. When fitting the dim phase spectrum with the NPEX plus Gaussian model, we find that the 
feature is better described in terms of an emission rather than an absorption profile. The statistical 
significance of this result, evaluated by means of an F-test, is between 2.91 x 10 -3 and 1.53 x 10~ 5 , 
taking into account the systematic errors in the background evaluation of HXD-PIN. We find that, 
the emission profile is more feasible than the absorption one for comparing the physical parameters in 
other phases. Therefore, we have possibly detected an emission line at the cyclotron resonance energy 
in the dim phase. 

Subject headings: X-rays: stars — X-rays: binaries — stars: pulsars: individual (4U 1626—67) - 
stars: magnetic fields 



1. INTRODUCTION 

Accreting binary pulsars are an ideal laboratory for 
studying the fundamental physics of radiative transfer of 
X-ray photons under strong (of the order of 10 12 G) mag- 
netic fields. A quite evident manifestation of such trans- 
fers are cyclotron resonance scattering features (CRSFs) 
observed in the spectra of many X-ray sources - — th e 
first being observed in Her X-l (iTru mper et al.lll978l) . 
CRSFs are caused by resonant scatterings between the 
Landau levels of electrons under the strong magnetic field 
near the surface of neutron star. The cyclotron reso- 
nance energy, E a , is strongly related to the magnetic 
field strength: 



S = 11.6S 12 -(l + ^)- 1 keV, 



(1) 



where B\i is the magnetic field strength in units of 10 12 
G, and z g represents the gravitational red-shift at the res- 
onance point. To date, CRSFs have been detected from 
16 accreting bi nary pulsars and measured their m agnetic 
field strength (|Miharalll995l: ICoburn et al.ll2002l) . How- 
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ever, the formation mechanism of CRSFs is not still fully 
understood due to the complex physics of resonant scat- 
tering and their continuum. Thus, the comparison of the 
CRSFs parameters observed in different sources, at dif- 
ferent luminosity and spectral states is important to gain 
information on the physical processes taking place. 

The Japanese X-ray satellite Suzak u, with its high 
sensi tivity hard X-ray detector (HXD; Takahash i et al.1 
[2001. has discovere d CRSFs from GX 304-1 
(|Yamamoto et al.l 120 111), success f ully d etected CRSFs 
from A0535+262 (jTerada et~"aLl 12006ft in its lowest 
luminosity epochs, an d performed detai led studies of 
the CRSFs of Her X-l (jEnoto et al.ll2008ft. 411 1907+0 9 
(jRlvers et al.ll2lTir]h and 1A 1118-61 rtSuchv et al]l20ll . 

4U 1626—67 was observed by Suzaku on 2006 March. 
The source is a 7.7 s pulsar in a low- mass X- 
ray binary system with an orbital period of 42 min 
(|Middleditch et alJH98ll : rChakrabartvill998l) . The X-ray 
sp ectrum shows low-energy emission lines from Ne and 
O (lAngelin i et al.lll995l: iSchulz et alJl2001t iKrauss et al.1 
I2007T ) and a feature at ~ 37 keV (jOrlandini et alJ l!998). 
The 37 keV absorption feature was discovered by Bep- 
poSAX and interpreted as the fundamental CRSF. The 
10-60 keV source flux is 4.4 x 10~ 10 erg cm~ 2 s _1 and 
the resonance energy corr esponds to a magnet ic field 
strength of - 3 x 10 12 G (jOrlandini et al.l [1991 . The 
broadband X-ray continuum exhibit s a strong phase de - 
penden ce as reported b y HE A O-l dPravdo et al.l|1979ft . 
Tenma (jKii et al.lll986h and RXTE ( |Coburnll2001l ). The 
cyclotron feature and the phase dependence variability 
provide information on radiative transfer in the plasma of 
the accretion column. In this paper, we used the March 
2006 Suzaku observation of 4U 1626—67 to investigate 
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the characteristics of the ~ 37 keV feature and the con- 
tinuum in the phase resolved spectra. The paper is or- 
ganized as follows: we describe the observations from 
Suzaku and the data reduction procedure in §2, show the 
results of the timing and spectral analyses in §3, and 
finally provide a discussion and summary in §4. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Suzaku observation of 4U 1626—67 

Suzaku is the fifth Japanese X-ray satellite 
(|Mitsuda et al.l 120071 ) to be launched and is equipped 
with two different instruments. One instrument is 
an X -ray imaging spectrometer (XIS; Kovam a et al.l 
120071) . a charge- coupled devi ce camera at the focu s 
of the X-ray telescope (XRT; ISerlemitsos et al.l 120071) . 
covering the 0.2-12 keV energy range. There are four 
XIS units; three of them (XIS-0, XIS-2 and XIS-3) 
are front-illuminated (FI) while the other one (XIS-1) 
is back- illuminated (BI) . The other instrume n t is a 
hard X-ray detector (HXD; iTakahashi et al.l 120071 ). 
which consists of PIN silicon diodes (HXD-PIN; 10-70 
keV) and Gd 2 Si0 5 Ce (GSO; 50-600 keV) scintillators. 
4U1626— 67 was observed by Suzaku on 2006 March 9 
UT 01:18 through March 11 19:38 at "XIS nominal" 
pointing position. The XIS operated in standard 
clocking mode and included the "1/8 window" option in 
order to have a time resolution of 1 s, without a charge 
injection function. 

2.2. Data reduction procedure 

We used data produced by Suzaku pipeline process- 
ing version 2.0.6.13 with the calibration versions hxd- 
20090515, xis-20090403 and xrt-20080709, and the soft- 
ware version HE AD AS 6.5.1. 

XIS and HXD events were screened by using standard 
criteria. We discarded events collected when the earth's 
limb was less than 5°, when the earth's day-night bound- 
ary angle from the line of sight was less than 25°, or when 
the spacecraft was in an orbit phase within 436 s (for 
the XIS) or 500 s (for the HXD) of the South Atlantic 
Anomaly ingress/egress. For XIS-1, we also excluded 
events in the time interval between 2006 March 10 UT 
06:58:30 and 07:20:54, because the cooling system was 
not properly working during this period. We accepted 
only XIS events with standard grades (0, 2, 3, 4 and 6) 
in the analysis. 

XIS spectra and light curves were extracted from a 
circle with a radius of 4.3 arcmin centered at the source. 
Spectra from the three FI XIS sensors were summed to- 
gether and the XIS intensity in the 0.4-10.0 keV energy 
band was ~ 5 c s" 1 for each sensor. We accumulated the 
XIS background spectra from a source-free region hav- 
ing the same exposure area as that of the corresponding 
sensors. The derived background spectra, exhibited an 
intensity of ~ 0.02 c s _1 for each sensor in the 0.4-10.0 
keV range, and is considered negligible. 

HXD screened events were used to obtain spectra and 
light curves for the PIN and GSO sensors. The HXD- 
PIN simulated non X-ray background (NXB) model 
was provided by the Suzaku HXD team, and we used 
the MET HOD = LCFITDT a nd METHODV = 2.0 
datasets (Fukaza wa et al.l [2009). The systematic un- 
certainty of the NXB model is estimated to be 2.31% 
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Fig. 1. — Energy-resolved light curves of 4U 1626—67 ob- 
served with Suzaku. Upper and middle panels show background- 
subtracted light curves obtained with XIS-0 in the 0.4-10.0 keV 
range and with the HXD in the 10—50 keV range, respectively. 
The bottom panel shows the hardness ratio of the above two pan- 
els. The gray solid line shows the average and the gray dotted lines 
show the ±3cr interval. 



for a 10 ks exposure in the 15-40 keV energy band 
(jFukazawa et al.ll2009D . Hereinafter, we assume 3% sys- 
tematic errors in the HXD-PIN analyses, and the cos- 
mic X-ray background (CXB) spect ra were sub tracted 
by using data acquired by HE A O-l (lBoldtlll987l) . How- 
ever, the source was not detected by the HXD-GSO. The 
NXB-subtracted HXD-GSO light curve and spectrum ex- 
ceeded the HXD-GSO NXB by ~ 1%, an d the latter were 
withi n typical systematic uncertainties (Fukazaw a et al.1 
I2009D . Net exposures were 102.7 ks for XIS-0, 2 and 3, 
101.4 ks for XIS-1 and 87.6 ks for the HXD. 

3. ANALYSIS AND RESULTS 
3.1. Timing analysis 

For the arrival time of each XIS and HXD-PIN event, 
barycent er correction was p erformed with the tool ae- 
barycen (jTerada et all 120081 ) . We surveyed the HXD- 
PIN data by using a folding technique to obtain a best 
spin period of -P(s U zaku) = 7.67795(9) s. This period is 
consistent with the spin-down trend that started to be 
observed from 1990 onwards, and that is calculated by 
using the spin-down rate of i> ~ — 4.8x 10 -13 Hz s _1 mea- 
sured with the Swift satellite over the period 2004-20 7 
(|Chakrabartv et al.lll997t ICamero-Arranz et alJl20Toh . 

XIS and HXD background-subtracted light curves of 
4U 1626—67 are shown in Figure 1. The light curves are 
plotted using a time bin of 7678 s; a multiple of P(s U zaku)- 
The XIS and HXD light curves have an intensity varia- 
tion of less than 7%, and the hardness ratio between them 
is approximately constant; the reduced chi-squared (x 2 ) 
value is about 1 when fitting a constant model to the 
hardness ratio. The average and their 3 a level are shown 
in the Figure 1. As a results, we confirmed that there 
is no significant spectral variation. Thus, X-ray spectra 
were extracted by using data over the entire observation. 



3 




Fig. 2. — Phase-averaged background-subtracted spectra of 4U 
1626—67 obtained with Suzaku XIS (black crosses denote XIS-FI 
data and light gray crosses denote XIS-BI data) and HXD (denoted 
by circles). Error bars are only for statistical predictions. The top 
panel shows the raw count rate energy spectra divided by those 
of Crab Nebula. The middle panel shows the spectra in counts/s 
keV - 1 along with the best-fit models (dotted lines) a blackbody, 
NPEX with Gaussian-absorption model and photoelectric absorp- 
tion (see text for details). Plots are presented without removal 
of instrumental responses. The bottom panel shows the residuals 
between the actual data and the best-fit model. 



3.2. Phase- averaged spectrum 

The top panel of Figure 2 shows the broadband 4U 
1626—67 phase-averaged raw count rate spectrum nor- 
malized by a Crab Nebula spectrum that was observed 
for calibration purposes. XIS source data were normal- 
ized by using a canonical model of the Crab Nebula spec- 
trum, since spectra from the Crab Nebula suffer from 
pile-up. HXD data were normalized by a HXD Crab 
Nebula spectrum observed on 2006 March 30 taken for 
the same observation parameters — position of the source 
with respect to the instrument FoV, grade criteria and 
PIN threshold file (ae_hxd_pinthr_20060727.fits)— as for 
4U 1626—67. The 4U 1626—67 source intensity is about 
20 mCrab at 20 keV, and the normalized spectrum is 
characterized by three main features: a soft X-ray ex- 
cess below ~ 3 keV , as reported by s e veral authors 
(lAngelini et all 119951: iSchulz et all 120011: iKrauss et all 
12007ft : a power-law continuum and a dip feature at ~35 
keV. 



Fig. 3. — The top panel shows the phase-averaged and 
background-subtracted spectra above 3.5 keV of 4U 1626—67 ob- 
tained with Suzaku XIS (FI: black crosses, BI: light gray crosses) 
and HXD (gray circles). The lower four panels show the residuals 
in units of standard deviations, Xi from fitting by: (a) the NPEX 
with GABS model without Fe-line emission model; (b) the NPEX 
with GABS model and Fe-line emission model; (c) the NPEX with 
CYAB and Fe-line emission model and (d) the ECUT with GABS 
model and Fe-line emission model. 



To model the phase-averaged spectrum, we fit the 
data with a ~ 0.24 keV blackbody for the soft ex- 
cess, a negative a nd positive power-law times exponentia l 
(NPEX) model (|Miharal U995t iMakishima eTail [19991) 
and a Gaussian- absorption (GABS) model (|Soong et al.l 
1990) for the continuum above 3 keV and a standard 
photoelectric absorption. The NPEX model 



f(E) = (A n E- a + A P E+P) exp ( - 



E 
kf 



(2) 



where E is the energy of an incident photon and kT is a 
plasma temperature (also called the cutoff energy), has 
two power law components with indices a and /3. The 
negative a term describes the spectrum at low energy, 
and the positive j3 term, which dominates as the energy 
increases, simulates the Wien peak in the unsaturated 
Comptonization regime when = 2.0. 

The best-fit model and its residuals are plotted in the 
lower panel of Figure 2. Above 3.5 keV only 1% of the 
NPEX negative term is contributed by the blackbody 
component. Therefore, hereinafter, we ignore XIS data 
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TABLE 1 

Best-fit parameters for phase-averaged spectra of 4U 1626- 



67 (d.o.f. = degrees of freedom) 



Parameter 


NPEX a 


NPEX with GABS a 


NPEX with GABS 


NPEX with CYAB 


ECUT with GABS 


a 


n 01 +0.03 
u - o± _0 03 


0.261°°* 


0.29i°°ol 


n o 9 +0.03 
u -°^_0 03 


0.651°;°? 


fcT(keV) 


r. o+O.l 


r.o±g:S 




7 2+ 1 - 1 




B c (keV) 










24.7l°'« 


E f (keV) 


— 


— 


— 







A n (xl0" 3 ) b 


o 1+0.4 




7 o + 0.3 

'• z -0.3 


7 o+O.l 
'•°-0.3 




A p (xl0- 5 ) b 


4 6+°' 3 
4 - D -0.3 


i fi +0.5 
1 -°-0.4 


i-8±8:l 






A(xl0~ 3 ) c 










10 0+ 139 
iu.u_ 75 


£ Fe (keV) 






6.62±°i° 


6.60l°;i; 


6. 60 (fix) 


a Fe (keV) 






o ? +0.13 


n c.q+0.13 


n o 7 +0.10 


EW Fe [eV] 






oo Q+16.1 

— 16.5 


or 9 +27.7 
od,z — 11.0 


31.3t 7 g 


E a (keV) 




fi +1.02 
°' - D -0.89 


o- 4 +l.O 


ok -+0.65 
a0 - '-0.55 


oo 4 +6.47 


OCRSF (keV) 




= 11+0.79 
°- 11 -0.68 


4-991°;^ 


5 56+ 1 - 77 
°- OD -1.08 


> 16.5 


T 




20 4+ 4 - 87 
zu ' 4 -4.20 


19.0if;° 


25.1±|;| 


7G 1+9.74 
'y- 1 _4.59 


X 2 1 d.o.f 


572.6/158 


182.1/155 


132.6/152 


134.6/152 


170.6/152 



a Without Fc-cmission line model. 

b Referring to equation (1), and defined at 1 keV in units of photons keV -1 cm -2 
c Normalization of the power-law. Defined at 1 keV in units of photons keV" 1 cm" 



below 3.5 keV and exclude the blackbody component 
from the fit. 

To evaluate the statistical significance of the CRSF in 
the phase averaged spectrum, first, we fit the spectrum 
with only NPEX continuum or NPEX x GABS model. 
The residuals of the NPEX with GABS model for en- 
ergy bands above 3.5 keV are shown in Figure 3(a), and 
the best-fit parameters are listed in Table 1. Second, 
we tes ted the significance using the F-test of iPress et al.l 
(|2007f ) routine. In this case, the F statistical value is 
defined as 

x?M (3) 



F = 



where the x\ an d x\ are statistic chi-squared, v\ and V2 
are degrees of freedom corresponding to the results of 
fittings using model 1 and model 2 (In this case model 
1 is only NPEX continuum and model 2 is NPEX with 
multiplicative component), respectively. The derived F 
statistic value is 3.1 which indicates that the probability 
of chance improvement of the x 2 is 3.7 x 10~ 12 . Thus, 
the CRSF is statistically significant feature in the phase 
averaged spectrum. We also found that the residuals 
in Figure 3(a) have a complex structure near the Fe- 
emission line energy around 6.6 keV. The fit improves by 
the addition of a Gaussian Fe-line model with an energy 
of 6.62 keV, a width (<Jp e ) of 0.37 keV and an equivalent 
width (EW) of 33.3 eV (Figure 3(b) and Table 1). These 
values are statistically consistent with the ASCA EW up- 
per limit of < 33 eV for an Fe-emis sion line in the 6.4- 
6.9 keV range ([Angelini et al.| [l995). The broad width 
may indicate the presence an Fe-line complex. If we use 
three narrow Gaussian Fe-lines, the Fe-line energies are 



6.37 



+0.07 



1+0.05 
-0.07 



and 6.92 



+ q q4 keV with correspond- 



0.08» 6 - 61 - 

ing EWs of < 11.1, < 12.2 and < 17.4 eV, although the 
fit is not improved significantly and we would not be dis- 
tinguish whether the Fe-emission line structure is made 
by three lines or not, due to the energy resolution of XIS 
(~150 eV). The X-ray unabsorbed fluxes in the 3.5-10 

keV and 10-60 keV energy bands are 1.27lg;5g x 10_1 ° 

l 

0.8 



erg cm 2 s 1 and 4.02 



+0.46 x 10 _10 erg cm -2 g - 



re- 



spectively, and the hard X -ray flux is consistent with 
the BeppoSAX observation (jOrlandini et al.l 11998). Ex- 
changing the GABS model with a cyclotron absorption 
model (CYAB; iMihara et all 19901 ) also fits the data well 
(Figure 3(c) and Table 1). 

If we replace the NPEX model by an exponen- 
tial cutoff power-law model (EC UT, which is the 
powerlaw*highecut model in Xspec: IWhite et aI1ll983| ) 
and fit the data with a GABS model, we obtain an ac- 
ceptable x 2 value (see Table 1), but substantial residual 
structures are seen above 20 keV (Figure 3(d)). Although 
the CRSF energy agr ees with BeppoSAX observations 
(jOrlandini et al.lll998f ). the best- fit continuum model is 
inconsistent with BeppoSAX results. 

3.3. Phase-resolved spectra 

Phase-res olved spectra acquired by HEAO-1 Tenma 
and R XTE (jPravdo et al.l[l979l : IKii et al J 11983 ICoburnl 
2001), show a strong phase dependence in the contin- 
uum emission of 4U 1626—67. Suzaku XIS and HXD 
folded light curves also indicate changes across the pulse 
(Figure 4). Specifically, while the 10-30 keV / 3.5-10 
keV ratio maintains the characteristics of the individ- 
ual light curves, the 30-50 keV / 10-30 keV ratio in- 
dicates spectral differences between phases 0.0-0.3 and 
the other phases. To investigate spectral changes across 
the pulse, we divided the events into eight phase-resolved 
spectra (Figure 5) and normalized these spectra by the 
Crab Nebula spectrum by using the method described in 
§3.2 (Figure 6). The Crab Ratios show that the ~ 20 keV 
peak is more prominent at higher pulse flux. Addition- 
ally, the Ratios flatten at lower pulse flux and the energy 
of resonance feature changes. The resonance features 
are particularly outstanding in the dim phase spectrum 
(</> = 0.125-0.250), where the hardness ratio (Figure 4) 
shows the large change in these features. 

We then performed spectral fitting to the phase- 
resolved spectra with the NPEX continuum and the Fe- 
emission line model in order to establish the pulse-phase 
dependence of the continuum and the CRSF. The Fe-line 
energy and the width of the Gaussian function were fixed 
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Phase 

Fig. 4. — Energy-resolved and background-inclusive pulse profiles of 4U 1626—67. Data were folded at P, 



(Suzaku) 



7.67795 s. Phase 



corresponds to MJD 53803.000092755. The average background count rates of 3.5-10.0 keV, 10.0-30.0 keV and 30.0-50.0 keV energy 
bands are 0.016 counts/s, 0.49 counts/s and 0.086 counts/s, respectively. 



TABLE 2 

NPEX CONTINUUM-ONLY BEST-FIT PARAMETERS FOR PHASE RESOLVED SPECTRA OF 4U 1626 — 67. 



pulse 


a 


kT 


A a 


4 a 
-rip 


EW Fe b 


X 2 / d.o.f 


phase 




(keV) 


(xlO- 3 ) 


(xl0~ 5 ) 


(eV) 




0.000-0.125 


-0 02+ ' 05 


4-4±C;| 


6 7 +0 ' 8 
'-0.7 


o q+2.3 
°-" i -1.9 


27+2 8 


81.2 / 81 


0.125-0.250 


u - 10 -o.n 


5.ii?;I 


t: O + 0.7 
b - B -0.5 


< 1.2 


8ll?I 


69.0 / 63 


0.125-0.250 


27+ - 02 
u - z '-0.02 


6 1+ 05 
D - L -0.4 


6 1+ 02 
D -- L -0.2 


(fixed) 




71.8 / 64 


0.250-0.375 


in+ 07 

u - iu -0.06 


4 q+0-2 


8.04 + o ? 




46+1 1 


93.4 / 97 


0.375-0.500 


U.01_q Qg 


5 0+ - 1 

°- u -o.i 




7 6+1-1 
'•"-0.7 


2ll« 


137.3 / 103 


0.500-0.625 


5fi +0 10 

u - OD -o.io 


= q + 0.1 

°-°-o.i 


10.4+1® 


7 9 +0.6 
'• z -0.6 


< 22 


297.9 / 129 


0.625-0.750 


64+ - 10 
u -°^-0.09 


5 2+ - 1 

°- z -o.i 




8.41,; 7 


35+ 5 


259.0 / 108 


0.750-0.875 


co+0.09 


= o + O.l 

-°-o.i 


9.811:5 




35±« 


247.2 / 91 


0.875-1.000 


4^+ - 07 


r. O+0.2 

-°-q.i 


9.9l ;° 8 


5.7 +°7 7 


oi +33 
01_ 27 


117.5/ 89 



Referring to equation (1), and defined at 1 keV in units of photons keV 1 cm 2 s . 

The energy and the width of the Fc-line arc fixed at 6.62 keV and (TFe — 0.37 keV, respectively. 



at values obtained from the averaged spectrum. The pa- 
rameters from the spectral fitting are listed in Table 2, 
and the ratios between the original data and the best-fit 
model are shown in Figure 7. We found that the ratio 
for the dim phase (cf> — 0.125-0.250) spectrum exhibits a 
monotone increase above 30 keV, and the residuals show 
evidence for the existence of an emission feature rather 
than an absorption (Figure 8(a)). In contrast, the other 
phases have a V-shaped ratio caused by a CRSF absorp- 
tion feature. Moreover, the fitting results for the dim 
phase spectrum show that the A p parameter is uncon- 
strained, strongly suggesting that the continuum spec- 



trum at the dim phase requires only the upper limit of the 
high-energy cutoff component. In other words, the result 
make the NPEX model effectively a power-law times ex- 
ponential cutoff model (cutoff pi in Xspec). The result 
is consistent with the trend that the 20 keV peak in the 
Crab Ratios is less prominent at lower pulse flux. We also 
tried other empirical models applied to accretion powered 
pulsar, ECUT and a power law with Fermi-Dirac cutoff 
(FDCO: lTanakalll986l) . and find that their best fit cutoff 
energies are consistent with 0. In other words, a power- 
law times an exponential cutoff model is the only one 
able to fit the dim phase spectrum. For this reason we 
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TABLE 3 

Best-fit parameters for phase resolved spectra of 4U 1626- 



■67, when fitting by NPEX with GABS model 



pulse 


a. 


kT 


A a 


A a 


E a 


a CRSF 


r 


EW Fe b x 2 1 d.o.f 


F/PCI 


C 


phase 




(keV) 


(xicr 3 ) 


(xlO- 5 ) 


(keV) 


(keV) 
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a Referring to equation (1), and denned at 1 keV in units of photons keV 1 cm 2 s . 

k The energy and the width of the Fc-line are fixed at 6.62 keV and <tf& — 0.37 keV, respectively. 

Referring to equation (3), and PCI is probability of chance improvement of the x 2 by adding a new component. 
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Fig. 5. — Phase-resolved and background-subtracted XIS (FI: 
black symbols below 10 keV, BI: gray symbols) and HXD (black 
symbols over 10 keV) spectra of 4U 1626—67 for phases </, = 0.000- 
0.125, 0.125-0.250, 0.250-0.375, 0.375-0.500, 0.500-0.625, 0.625- 
0.750, 0.750-0.875 and 0.875-1.000 from top to bottom, offset by 
multiples of 10 — 2 . Arrows denote the upper limit. 
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Fig. 6. — Ratios between the phase-resolved and background- 
subtracted XIS and HXD spectra of 4U 1626—67 and those of the 
Crab Nebula. The XIS source data were normalized using the 
canonical model of the Crab Nebula since spectra from the Crab 
Nebula suffer from pile-up. Phases are indicated in the same man- 
ner as Figure 5, but are offset by multiples of 10 . Arrows denote 
the upper limit. 



fit the dim phase spectrum by using the NPEX model 
with fixed A p = 0. As we can see in Table 2, the differ- 
ence between the x 2 values for freely selected A p and for 
fixed A p = is negligible. Therefore, we henceforth fix 
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Fig. 7. — Ratios between the data and best-fit NPEX-only model 
for each phase-resolved spectra. Phases are indicated in the same 
manner as Figure 5. Arrows denote the upper limit. 
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Fig. 8. — Upper panel shows spectra in the dim phase (<j> = 0.125— 
0.250) of 4U 1626-67 obtained with Suzaku XIS (FI: black crosses, 
BI: light gray crosses) and HXD (circles). Triangles denote the 
reproducibility of NXB for HXD-PIN. Histograms in the upper 
panel are the best- fit by NPEX continuum (A p = 0). The lower 
three panels show residuals, \, when fitting by (a) NPEX (A p = 0) 
continuum-only; (b) NPEX (Ap = 0) with GABS model and (c) 
NPEX (A p = 0) and Gaussian models. 



A p = for the NPEX model when fitting the dim phase 
spectrum. 

Next, we included the GABS model and refit the phase- 
resolved spectra. The resulting parametrized models fit 
all of the phase-resolved spectra well (see Table 3 and 
Figure 9), and residuals above 30 keV in the dim phase 
were eliminated (see Figure 8(b)). We also evaluate the 
significance of including a GABS component using the 
same method discussed in section 3.1. The derived prob- 
abilities are summarized in Table 3, that shows that the 
absorption factor is not statistically significant in the 
dim phase. This significance did not change even when 
we adopted the CYAB model (as opposed to the GABS 
model as the cyclotron feature model) for spectral fitting. 

Since the GABS model reproduces an absorption fea- 
ture, we fit the phase-resolved spectra by the NPEX plus 
a Gaussian model in which we allowed the Gaussian nor- 
malization to be both positive and negative. The de- 
rived parameters are shown at Table 4 and Figure 10. 
The results indicates that the emission line produces the 
residuals from the NPEX fitting in Figure 8(a). We can 
evaluate how much the additional ter m improved the fi t 
using the F-test for additional terms (|Bevingtonl 11969). 



The F statistical value, hereafter called F x , is defined as 



(4) 



Each suffix denotes model 1 and 2, NPEX only and 
NPEX + Gaussian model, respectively. The derived 
probability from F x also summarized at Table 5. The 
results show that the data support the presence of the 
emission line in the dim phase with F x of 7.688, cor- 
responding to a probability of chance improvement of 
the x 2 °f 1-9 x 10~ 4 . The residuals of the fitting are 
shown in Figure 8(c). The obtained resonance energy, 
E a = 40.5l 1 47 keV, is statistically consistent with those 
observed in the other phases. 

To examine the chance probability of improvement 
of the derived F- value, we simulate th e spectra follow- 
ing a procedure similar to section 5.2 of iProtassov et all 
(2002). First, we simulate 10,000 data sets using "fakeit" 
in XSPEC Versionl2 according to the NPEX (A p = 0) 
model with parameters fixed at best fit values listed in 
Table 2, same instrumental response, background files 
and same exposure time. Each data set is binned us- 
ing same binning as that used to real data. Second, we 
fit the simulate data using only NPEX (A p = 0) and 
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Fig. 9. — Fit parameters from as a function of pulse phase for 4U 1626 
plotted. 
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Best-fit values by NPEX with GABS model (Table 3) are 



TABLE 4 

Best-fit parameters for phase resolved spectra of 4U 1626—67, when fitting by NPEX plus Gaussian model. 
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a Referring to equation (1), and denned at 1 keV in units of photons keV cm s 
k Normalization of the Gaussian is defined in units of photons cm -2 s _1 . 

c The energy and the width of the Fe-line are fixed at 6.62 keV and <7p e — 0.37 keV, respectively. 

Referring to equation (4), and PCI is probability of chance improvement of the % 2 by adding a new component. 
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Fig. 10. — Fit parameters from as a function of pulse phase for 4U 1626—67. Best-fit values by NPEX plus Gaussian model (Table 4) 
are plotted. 
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Fig. 11. — Observed value of the F x plotted as dashed line against 
probability histogram of F x obtained by 10,000 simulations (See 
the text). 

NPEX = 0) + Gaussian model allowing their nor- 
malization to be negative. Third, we calculate the F x 
derived by these two model fittings. Last, we evaluated 
the number of simulated spectra with F x exceeding the 
7.688, which is derived by fitting of real data to estimate 
the approximate probability. The derived histogram of 
the F x is shown at Figure 11 and the evaluated chance 
probability is 3.0 x 10 -4 which is consistent with the F x 
value obtained from real data using equation (4). In or- 
der to evaluate the effect of systematic error of HXD-PIN 
non-X-ray-background model, we also refit the dim phase 
spectrum taking account into the 3 % systematic errors 
of HXD-PIN. The F x is 7.688±|;||g and the probability 
is between 2.91 x 10 _3 and 1.53 x 10~ 5 . Therefore, even 
though the NXB level of HXD-PIN is underestimated, we 
have possibly detected the emission line at the cyclotron 
resonance energy in the dim phase. Although there is a 
~ 85 % possibility that the structure in the dim phase is 
produced by absorption feature, we stress that the con- 
tinuum requires a hard component with a large cutoff 
energy, in addition to a wide and deep absorption fea- 
ture, as seen in Figure 12. Furthermore, the kT value 
at the dim phase matches those at the other phases with 
the emission model, whereas a large jump in the value is 
required with the absorption model, as shown in Figure 
9. 

4. DISCUSSION 

We observed the 4U 1626—67 pulsar with the Suzaku 
satellite in March 2006. The pulse period deter- 
mined by the HXD-PIN of Suzaku is P(s U zaku) = 
7.6 7795(9) s, and this va l ue ag rees with that reported 
by iCamero-Arranz et al.l (|2010f) . In the pulse phase- 
averaged spectra, the fundamental CRSF is clearly de- 
tected at 37.4to g keV (GABS model) and the contin- 
uum emission is better reproduced by the NPEX model 
than the ECUT model. Phase-resolved spectral analy- 
ses consistently show strong phase dependence, both in 
the continuum emission and in the CRSF, particularly 
in the dim phase. Quantitatively, the positive compo- 
nent of the NPEX model vanishes near the dim phase, 



and the CRSF in the dim phase is reproduced by using 
a Gaussian shaped emission line model. T he statisti- 
cal sig nificance of the emission line using the iBevingtonl 
(1969) F-test routine is 2.91 x 10~ 3 and 1.53 x 10" 5 tak- 
ing account into the systematic errors in the background 
evaluation of HXD-PIN. The results means that even if 
the the real NXB level of HXD-PIN is 3% higher, which 
is the worst case of underestimating NXB level, the sig- 
nificance level is still around 3a level. Thus, we have 
possibly detected the emission line at the cyclotron res- 
onance energy in the dim phase. 

In this section, we discuss the nature of the accretion 
column on the magnetic pole of the neutron star that pro- 
duces a cyclotron emission line at the dim phase rather 
than an absorption feature. 

What is the o rigin of the emission feature in the dim 
phase? iNagej (|1981[ ) pointed out that if the photons 
originate from an optically thick accretion column (with 
a Thomson optical depth tt > 1) without temperature 
and density gradients, the CRSF is expected to be an ab- 
sorption feature with Comptonized spectra. In contrast, 
if the photons originate from an optically thin accretion 
column (tt < 1), the CRSF is expected to be an emis- 
sion feature with a free-free continuum. To compare this 
theoretical model with the Suzaku spectra, the unfolded 
spectra in the on-phase (0 = 0.625-0.750) and the dim 
phase ((f) — 0.125-0.250) are shown in Figure 13. Quali- 
tatively, our observed results are consistent with the the- 
oretical predictions; the Wien hump disappears with the 
appearance of an emission-like feature. This resemblance 
suggests that we observed the cyclotron resonance scat- 
tering emission feature caused by collisional excitation of 
a Landau level, followed by radiative deexcitation. The 
deciding factor as to whether we observe absorption spec- 
tra or emission spectra is thus the optical depth of the 
accretion column, and for the case of 4U 1626—67, since 
change occurs according to the spin phase, changes in 
optical depth arise from changes in viewing angle. 

The optical depth, r, is generally expressed as 

t = n e al. (5) 

Here, n e is the electron number density, a is the scat- 
tering cross section and / is the plasma length-scale 
through which a photon passes. In a strong magnetic 
field, where the strength of the field is greater than the 
thermal energy, the scattering cross section depends on 
the direction of pho tons relative t o the magnetic field 
vector. According to iHeroldl (|1979f) . the scattering cross 
section of photons averaged over the energy and the (or- 
dinary and extraordinary) modes is described by 

when the photons have energies far less than E a , and they 
propagate parallel to the magnetic field. Here, ot is the 
Thomson scattering cross section and E is the mean pho- 
ton energy. Thus, the optical depth propagating parallel 
to the magnetic field, n , is 

In contrast, the averaged cross section for photons prop- 
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Fig. 12.— The dim phase vFv spectra of the 4U 1626-67 obtained with XIS and HXD. The NPEX-only (A p 
the histogram. (Left) When fitted by NPEX (A p = 0) with GABS model and (right) when fitted by NPEX (A p = 



= 0) models are shown in 
0) and Gaussian models. 
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Fig. 13.— The unfolded spectrum of 4U 1626-67. Crosses de- 
note on-phase spectra (</> = 0.625-0.750). Squares denote off-phase 
spectra {<j> = 0.125-0.250). 



agating perpendicular to the magnetic field, a±, can be 
approximated by 

<t±_ w a T , (8) 

and the optical depth propagating perpendicular to the 
magnetic field, tj_, is 



(9) 



Therefore, from equations (7) and (9) we see that the op- 
tical depth is anisotropic, and the probability of Comp- 
tonization of photons depends on the direction that 
the photons are propagating through the plasma. This 
anisotropy of optical depth in the Comptonization pro- 
cess is expected to cause the phase dependence of the 
Wien hump. 

To verify the phase dependence of the Wien hump, we 
again fit the phase-resolved spectra by the NPEX model, 
except for the dim phase, where the kT value are instead 
fixed at the phase-averaged value (Table 1) by assuming 
thermal equilibrium. Here, the CRSF is described by 
the GABS model and the spectral fits are found to be 
acceptable. Figure 14 shows the background-subtracted 
pulse profile for the 10-30 keV energy band and compares 
this with values of A p plotted as a function of the pulse 
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Fig. 14. — The plotted line is the background-subtracted pulse 
profile in the 10-30 keV energy band. The circles denote the best-fit 
values of the normalization of the positive power-law term (i.e., the 
Wien hump) when fitting by NPEX model, presented as a function 
of the pulse phase (see text for details). Values are normalized by 
each maximum value. 



phase. These two plots clearly agree. In other words, X- 
ray fluxes in the 10-30 keV range are represented by only 
the Wien hump. Since, as described above, radiation is 
caused by the optical depth of scattering, the dim phase 
corresponds to the optically thin phase. Therefore, this 
result is explained correctly by our prediction. 

To estimate the value of rn qualitatively, we use the 
simple assumptions that the plasma is at thermal equi- 
librium and that the accretion column is homogeneous. 
Under these conditions, E is equal to the electron tem- 
perature. Therefore, the ratio between T\\ and t± is es- 
timated by using the observed kT ~ 6.71 keV and the 
E a ~ 37.4 keV derived from fitting the phase-averaged 
spectra by the NPEX with GABS model (Table 1): 



T|| 



i± 



kT\ 2 l\\ in 

_)_I = 3 .2xlO-*JL. (10) 
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Thus, if the difference between Z|| and l± is small (Zii ~ 
the optical depth for the dim phase spectra is ap- 
proximately one or two orders of magnitude lower than 
that in the bright phase. If we assume that the Thom- 
son optical depth of the bright phase is around 10, the 
cyclotron resonance scattering absorption feature can be 
seen as an emission feature at a certain viewing angle. 
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theoretical studies that take into account emission power. 
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